Atomic current across an optical lattice 
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We devise a microscopic model for the emergence of a collision-induced, fermionic atomic cur- 
rent across a tilted optical lattice. Tuning the - experimentally controllable - parameters of the 
microscopic dynamics allows to switch from Ohmic to negative differential conductance. 
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Ultracold atoms in optical potentials provide a versa- 
tile tool for the experimental study of many-body dynam- 
ics, with an unprecedented degree of accuracy Q, 0, 
Not only do these systems allow the faithful experimental 
realization of fundamental model Hamiltonians of solid 
state theory, and, thus, lend themselves as efficient ana- 
log quantum simulators Beyond that, they bridge the 
gap between the single particle and the thermodynamic 
limit: if not now, so certainly in the near future experi- 
mentalists will be able to load a precisely controlled num- 
ber of particles into engineered potentials, with particle 
numbers tunable from one to intermediate or large val- 
ues 0. Hence, by continuously varying the system size, 
we will be able to control the emergence of macroscopic, 
thermodynamic observables from perfectly deterministic 
Hamiltonian dynamics, in experiment and theory. This 
holds the potential of a much deeper and quantitative un- 
derstanding of quantum statistical laws, together with a 
variety of possible applications in the context of the ever 
advancing miniaturization of modern technology @ . 

In the present contribution, we address a specific prob- 
lem of exemplary importance within this context: the 
emergence of a macroscopic current across a periodic po- 
tential under static forcing. It is well known that elec- 
trons in a perfectly periodic lattice exhibit coherent Bloch 
oscillations, under static forcing, and that no net trans- 
port, i.e., no current across such lattice occurs in the 
absence of any (incoherent or dephasing) relaxation pro- 
cess ■ Only the latter will induce a net drift velocity of 
the electrons, giving rise to a measurable current. While 
relaxation processes due to phonon scattering or impu- 
rities are abundant yet largely uncontrolled in solids, a 
clean and perfectly controllable realization of this funda- 
mental transport problem is feasible in a quantum opti- 
cal setting, where we can engineer the environment: By 
loading spin polarized fermions together with ultracold 
bosons into a one dimensional optical lattice, we estab- 
lish a strictly analogous, idealized scenario. The bosons 
act as a bath for the fermions, with relaxation induced 
by fermion-boson collisions. 

Our results are obtained from a microscopic model 
for very moderate system sizes - typically less than ten 
bosons interacting with a single fermion, on approx. ten 
lattice sites. Yet, due to exponential increase of the 
many particle system's Hilbert space dimension with bo- 



son number N and lattice size L, this is sufficient to enter 
the regime of asymptotic convergence into the thermo- 
dynamic limit Thus, we describe how a macroscopic 
current emerges from perfectly Hamiltonian microscopic 
dynamics, with all ingredients under essentially perfect 
experimental control. 

Our model Hamiltonian reads 

H = H F + H B + H int , (1) 

and decomposes into the (single particle) fermionic part 

H * = ~t (E \ l + + h - c ^j + ™E iwi , (2) 

the (many particle) bosonic part 

H B = + h - c -"j + ~Y E fll( - 7il ~ ^ ' 

(3) 

and a term which mediates the collisional interaction be- 
tween fermions and bosons, 

H int =W FS ^2h l \l){l\. (4) 
i 

Since we will be interested in the fermionic dynamics, 
the fermionic tunnelling coupling Jp in (2) sets the en- 
ergy scale, and, correspondingly, h/ Jp is the natural unit 
of time. Equations (1-4) describe the dynamics of ultra- 
cold atoms in the lowest band of a ID lattice, assuming 
short range (on-site) interaction, and tunnelling coupling 
only between adjacent lattice sites I. Low temperatures 
as well as moderate field strengths Fd < Jf are required 
for consistency with the single band approximation. The 
single particle ansatz for the fermionic part is justified 
since the Pauli principle forbids occupation of the same 
site. Our model additionally implies that only fermions 
experience a static forcing F, while the bosons see a pe- 
riodic potential (period d) without static forcing. Such 
situation can be realized by a suitable choice of the inter- 
nal electronic states the fermions and bosons are initially 
prepared in |9(. 

We impose Jb = Jf for the bosonic component and 
focus on filling factors n — N/L ~ 1 with N bosons dis- 
tributed over L lattice sites (for more than one fermion, 
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N represents the boson number per fermion). These fill- 
ing factors interpolate between single-particle dynamics 
and the thermodynamic limit, for the lattice sizes we con- 
sider. The above choice of parameters is in reach for state 
of the art experiments: Jb — Jf/2 has already been re- 
alized in experiments 0] on a mixture of 40 K and 87 Rb, 
with Jp ~ 0.06-Er, and experimental techniques for the 
relative adjustment of Jp and Jb are available P.I^.Hl|. 
Here, En = h 2 k 2 /2m is the single photon recoil energy 
seen by the 40 K atoms. With the above value of Jp, 
the relevant time scale H/Jp ~ 0.3 s is easily resolved. 
Finally, for our subsequent theoretical/numerical treat- 
ment, the total dimension J\f of the Hilbert space scales 
as TV = LJ\f B , where M B = (N + L - l)\/N\{L - 1)! is 
the dimension of the bosonic subspace. The exponential 
increase of J\f with Nb and L is thus evident. 

In the absence of interactions between the fermionic 
and bosonic component (Wfb = in Q), the fermions 
will undergo coherent Bloch oscillations with a period 
Tb = 2-Kh/Fd 0| . However, finite values of Wfb will 
induce a collision-induced relaxation of the fermionic dy- 
namics, eventually resulting in a nonvanishing current 
across the lattice. For weak coupling between fermions 
and bosons, we can describe the reduced fermionic dy- 
namics in second order perturbation theory in Wfb • Fur- 
ther we use the fact that the bosonic subsystem is known 
to exhibit chaotic level statistics (in the random matrix 
sense), over a broad range of Wb and Jb [S 0, and 
therefore can substitute for a noisy environment |13| . 
provided its initial state is characterized by equally dis- 
tributed populations of the different sublevels (this is, in 
the high-temperature limit). Indeed, the boson- number 
cross-correlation function 



Rl,m(t,t') = Tr[hi(t)h m (t')] , 
fii(t) = exp(iHBt/h)fiiexp(—iHBt/h) 



(5) 



is very well approximated by the cross correlation gen- 
erated by a random matrix Hamiltonian (chosen from 
the Gaussian Orthogonal Ensemble 0]), as illustrated 
in Fig. 1 for (i,0) and R h2 (t,0). For TV > 2 and L 
such that Nb > 100, we numerically extract a well de- 
fined, finite time scale r ~ 3H/Jb for the decay of these 
correlations |l7|. 

Consequently, if we choose lattice constant d and static 
field strength F such that r <C Tb 0], the correlation 
function of the bath can be approximated by a ^-function 



Rl,m(t,t') = n 2 Si. m S 



t-t' 



(6) 



on the typical time scale of the fermionic dynamics, given 
by the Bloch period Tb- This justifies a Markov approx- 
imation after tracing out the bosonic degrees of freedom, 
and we end up with the master equation for the fermionic 
density matrix, 



dt 



H F (t),p 



{F) 



l.m 



7(1 - 5l, m )p. 



{F) 



(7) 
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FIG. 1: Time evolution of the number correlation functions 
Ri,i(t, 0) (top) and Ri^(t, 0) (bottom) of the bosonic bath 
©, initially prepared in an equally weigthed superposition of 
the bosonic eigenstates. N = 7, L = 9. The ratio between in- 
teratomic interaction strength and tunnelling coupling is fixed 
at Wb/Jb = 3/7, what insures that we are in the parame- 
ter range of Wigner-Dyson spectral statistics, at n ~ 1 @. 
Time is measured in units of the site-to-site tunnelling period 
tb = hi Jb- The dashed lines show the bosonic correlation 
functions which result from substituting the Bose-Hubbard 
Hamiltonian © by a random matrix of the Gaussian Or- 
thogonal Ensemble. 



where the decay rate is given by 



rn 2 W 2 B 

h 2 



3n 2 W 2 B 

hJ B 



(8) 



An analytic solution [T^ of (0) leads to the following 
expression for the mean fermionic velocity in the lattice: 



v(t) = v e 7 sin(w^) , vq 



Jpd 



(9) 



in perfect agreement with the numerical result obtained 
from a direct propagation of the dynamics, with the 
bosonic bath prepared in an equally weighted superpo- 
sition of the eigenstates of the Bose-Hubbard Hamilto- 
nian, garnished by random phases (tantamount of the 
high temperature limit), while the fermionic subsystem 
is launched in a Bloch wave with vanishing quasimomen- 
tum. As illustrated in Fig. 2, the fermionic Bloch oscil- 
lations decay irreversibly, with a decay constant 7 which 
- in agreement with Eq. © - is fully controlled by the 
filling factor n, by the fermion-boson interaction strength 
Wfb, and by the bosonic tunnelling coupling Jb- 

However, Fig. 2 indicates no nonvanishing average drift 
of the fermions across the lattice, hence no current is ob- 
served. This is a consequence of our above derivation: 
while the master equation description provides a quite 
satisfactory and quantitative picture of the collision in- 
duced decoherence of the fermionic Bloch oscillations, it 
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FIG. 2: Normalized mean velocity of the Fermi atoms, in units 
of no = JfJ/Fi, for Fd = 0.57 x Jf, Jb = Jf, for fermion- 
boson interaction strengths IFfb = 0.101 x Jf, 0.143 x Jf, 
0.202 x J F (from top to bottom), and N = 7, L = 9. The typi- 
cal time scale of the interaction induced decay of the fermionic 
Bloch oscillations fits the time scale 7 -1 predicted by eq. (JHJ 
(dash-dotted vertical lines) very well. The same observation 
holds for fixed Wfb and variable Fd (not shown here). 



also implies that there is no back-action of the fermion- 
boson coupling on the bath (this is the essence of the 
Markov approximation). Consequently, no net energy 
can be transfered from the fermionic into the bosonic 
component (the latter being initially in a fully thermal- 
ized state), as clearly manifest in the time evolution of the 
average energy of the bosonic component, represented by 
the dashed line in the bottom panel of Fig. 3. If, instead, 
we choose a bosonic initial condition given by a low tem- 
perature Boltzmann distribution P(ej) ~ exp(— ej/fcsT), 
with temperature ksT ~ 2.86 x Jb uM, over its energy 
levels ti, the bath can be heated by the collisions with 
the fermions, and thus extract energy from the fermionic 
component. Such choice of the bosonic initial state leads 
to the relaxation of the fermionic Bloch oscillations, as- 
sociated with a (non-Markovian) energy increase of the 
bosonic component, see Fig. 3. Correspondingly, the 
fermionic mean velocity exhibits a nonvanishing drift, on 
top of the oscillatory behaviour due to the Bloch dy- 
namics. A finite fermionic current across the lattice is 
observed in the upper panel of Fig. 3! Since we will be 
dealing with a closed system of finite size (finite particle 
number and finite lattice length), the bosonic bath has a 
finite heat capacity and cannot act as a reservoir over ar- 
bitrary time scales. Consequently, the relaxation-induced 
current ceases as soon as the bath is fully thermalized by 
draining energy from the fermionic dynamics. 
Fermionic drift 

v(t) = v(t) - v e~ 7 * sin(cj B t) (10) 

and bosonic energy gain AEb are related through the 



FIG. 3: Mean velocity v(t) of the fermionic component (top, 
solid line), in units of vo = Jpd/h, and average energy Es(t) 
of the bosonic component (bottom, solid line), in units of 
Jf, for a low temperature (fesT ~ 2.86 x Jb) initial state 
of the bosonic bath. Fd = 0.143 x J F , IFfb = 0.143 x J F , 
Wb/Jb = 3/7, J B = Jf, N = 7, L = 9. The dashed lines 
indicate the corresponding solutions for a thermalized initial 
state of the bath, ksT ~ 150 x Jb- For a nonequilibrium 
initial condition of the bath the finite drift velocity v(t) man- 
ifests as a clear, nonvanishing offset of the fermionic mean 
velocity with respect to the Markovian result indicated by 
the dashed line (top). The dash-dotted line in the bottom 
plot is obtained by integration of eq. (II IB . 



classical relation AEb = FAx = FvAt, i.e., 



dE E 
dt 



Fv{t) 



(11) 



Integration of Ijllfl with v as defined in 110|) , together with 
v(t) from the exact solution of the Schrodinger equation 
under QJ, leads to a very good fit of the exact time de- 
pendence of the bosonic mean energy, in the lower panel 
of Fig. 3. We can thus deduce the current-voltage charac- 
teristics for the fermions, i.e., the dependence of the drift 
velocity v (equivalent to the current, modulo the carrier 
density) on the static force (or potential difference across 
the lattice) which generates the current: For sufficiently 
short times, eq. Ijllf) can be rewritten as AEb = FvAt, 
where Fv can be extracted from the exact result for 
Es(t) as displayed in Fig. 3 (averaging over the resid- 
ual Bloch oscillations during the first two Bloch cycles). 
Dividing this initial growth rate of Es(t) by F, we obtain 
the dependence of the "current" v on the "voltage" F in 
Fig. 4, which exhibits a marked transition from Ohmic 
(v ~ F) behaviour to negative differential conductivity 
(v ~ 1/F), as F is increased. Remarkably, this faith- 
fully reproduces the qualitative behaviour predicted by 
the phenomenalogical Esaki-Tsu relation [ljj for a cur- 
rent across doped superlattices, 
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FIG. 4: Drift velocity v (stars), in units of Vo = Jpd/h, vs. 
static forcing Fd, in units of Jf. The Bloch frequency ujb ~ 
Fd increases with increasing F, while the effective decay rate 
7 ~ 0.037 J F /h, derived from © for TV = 7, L = 9, IF FB = 
0.143 x Jf, Wb/Jb = 3/7, Jb = Jf, is constant. Note the 
clear transition from Ohmic behaviour (v ~ Fd) to negative 
differential conductance (v ~ 1/Fd), at a finite value of Fd. 

with quantitative differences, notably for large F, due 
to the finite size of our atomic sample. Note that, in 
contrast to experiments in semiconductor superlattices, 
the effective scattering rate 7 is perfectly controlled and 
continuously tunable (e.g., via Feshbach resonances 0|) 
in our present quantum optical setting, as a consequence 
of eq. ©• 



In conclusion, we have shown that a directed current 
of spinpolarized ultracold fermionic atoms across a one 
dimensional optical lattice under static forcing can be 
induced by collisional interaction with a bosonic admix- 
ture. While the current increases linearly with the static 
force in the limit where the collision induced decay rate 
(independent of the static forcing!) is much larger than 
the fermionic Bloch frequency (increasing linearly with 
the static field strength), it decreases inversely propor- 
tional to F in the opposite limit. This crossover has an 
intuitive cause: Finite, collision induced decay rates in- 
duce diffusive transport at arbitrarily weak forcing (many 
collision events during one Bloch cycle), while strong 
forcing and small collision rates (few collisions during 
many Bloch cycles) essentially reestablish the Bloch os- 
cillations, and suppress the directed current. 

Our theory, based on a microscopic, strictly Hamilto- 
nian picture of the many-body dynamics, is amenable to 
experimental scrutiny. It also identifies the microscopic 
origin of first experimental observations which already 
pointed in this direction, on the basis of a purely phe- 
nomenological interpretation Finally, our model al- 
lows for bosonic baths composed of an arbitrary number 
of particles. Future work will explore the limit of very 
small baths, with small Hilbert space dimensions and a 
reduced density of states. In this regime, non-Markovian 
effects due to the granularity of the bath are expected, 
which will be in reach for direct experimental probing. 
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